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FOREWORD 


This  report  was  prepared  bjr  Leon  J,  Bowser  as  au  In-house  effort 
under  Projeot  Number  1367,  'Struotural  Deilgpii  Orittria*t  Tiosk  Number 
136702.  "lerospaoe  Vehicle  Structural  Design  Oriterla'*  Iho  vcrk  vad 
aoooQplished  under  the  direction  of  C«  J.  Sohmldi  project  engineer* 

This  report  covers  work  performed  from  February  I961  throo£^  DfO* 
eober  1961. 
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the  report  includes  a  detailed  account  of  some  Aeronautical  Systems 
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duction  to  spectral  analysis  techniques. 
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Introduction 


This  rejxjrt  contains  a  review  and  analysis  of  the  state-of-the-art 
status  regard inc5  application  of  power  spectral  density  techniques.  This 
report  does  not  presume  to  supersede  any  of  the  literature  in  the  areas 
of  generalized  harmonic  analysis,  transform  calculus,  or  randan-process 
theory.  Stated  differently,  this  report  is  concerned  solely  with  the 
application  and  exploitation  of  established  methods  and  techniques,  and 
not  the  definition  or  exposition  of  fundamental  concepts  of  post-doctor¬ 
ate  level  mathematics.  It  is  fully  acknowledged  that  a  philosophical 
analysis  of  these  techniques  is  a  legitimate  study.  However,  In  applica¬ 
tions,  the  abstract  mathematical  models  serve  as  tools  and  different 
models  can  describe  the  same  empirical  situation.  The  manner  in  which 
mathematical  theories  are  applied  does  not  depend  on  the  preconceived 
ideas;  it  is  a  purposeful  technique  depending  on,  and  changing  with 
experience . 


A  I“toral 

■To  analyze  time  series  effectively  we  must  do  the  same  as  in  any 
other  area  of  statistical  technique  “Fear  the  Lord  and  Suuna  the  Devil" 
by  admitting  that: 

(1)  The  complexity  of  the  situation  we  study  is  greater  than  the 
complexity  of  that  description  of  it  offered  by  our  estimates. 

(2)  Balancing  of  one  ill  against  another  in  choosing  the  way  data 

is  either  to  be  gathered  or  to  be  initially  analyzed  always  requires  know¬ 
ledge  of  quantities  which  cannot  be  merely  hypothesized  and  which,  in  many 
cases,  we  cannot  usefully  estimate  from  a  single  body  of  data,  such  as 
ratios  of  (detailed)  variance  components  or  extents  of  non-nomality . 
Theoretical  optimizations  based  upon  specific  values  of  such  quantities 
may  be  useful  guides.  However,  theoretical  optimizations  should  be  used 
only  when  the  failure  of  past  experiences  and  the  present  data  to  give 
precise  values  for  theso  quantities  is  recognized  and  allowed  for. 

(3)  There  is  no  substitute  for  repetition  as  a  basis  of  assessing 
stability  of  estimates  and  establishing  confidence  limits. 

(4)  Asymptotic  theory  must  be  a  tcjl,  and  not  a  master. 

The  only  difference  is  that  one  must  be  far  more  conscious  of 
these  acceptances  in  time  series  analysis  than  in  most  other  statistical 
areas."  Some  concrete  examples  of  these  principles  are  briefly  discussed 
in  this  section  and  expanded  Ufon  in  subsequent  sections  and  appendices. 


♦Dr  John  W,  Tukey,  Princeton  U.  Tech  Rpt  No.  91,  I«3arch  I961. 


In  rcany  structural  design  problens,  the  knov/lodge  of  the  physical 
phenomena  is  not  precise  enough  to  justify'  exact  predictions  with  respect 
to  each  individual  obsei’vation.  On  the  ons  hand,  the  loading  is  never 
exactly  known,  and  on  the  other,  the  physical  system  may  bo  so  cccipli- 
cated  that  detailed  calculations  for  the  prediction  of  dynamic  reaponsei 
are  difficult,  (e.g.  The  random  nature  of  the  atmospheric  turbulence 
as  shown  in  the  records  is  duo  to  the  influences  of  viscosity,  pcressure, 
density,  temperature,  humidity,  and  velocity  distributions,  varying  in 
such  a  complicated  manner  that  mechanistic  predictions  based  on  the  hydro- 
dynamic  equations  are  practically  impossible.) 

In  such  cases  the  physical  data  can  be  regarded  as  a  sot  of  statisti¬ 
cal  data  of  random  processes.  Logical  statements  of  the  forcing  functions 
are  therefore  otati3tical--not  mathematical  but  statistical.  This  is 
very  important  here  because  a  mathematical  statement  usually  makes  an 
assertion  which  is  true  in  each  individual  case,  and  if  there  is  one  case 
which  violates  the  mathematical  statement  it  also  disproves  the  same 
etatemsnt.  V/hilo  a  statistical  statement  ensures  only  that  over  the  long 
run  its  truth  will  provide  itself.  For  example,  a  typical  statistical 
statement  may  be  that  the  probability  of  obtaining  heads  in  a  single  toss 
Of  a  coin  is  1/2.  However,  two  consecutive  tails  surely  is  not  sufficient 
to  destroy  faith  in  the  statistical  hypothesis. 

If  one  recognizes  then  the  uncertainty  of  predicting  the  forcing 
functions  in  a  given  problem,  the  problem  immediately  becomes  statistical, 
and  one  inuot  attempt  to  state  the  main  features  of  the  responses,  as  well 
as  those  of  the  forcing  functions  in  terms  of  statistical  averages  and 
probability  distributions.  Another  important  example  is  structural  fatigu® 
under  random  repeated  loading.  Introductory  material  relevant  to  this 
problem  is  contained  in  Appendix  II  of  this  report.  • 


An  excellent  sunaation  of  the  state-of-the-art  aa  of  I-lay  1957  la 
afforded  by  Dr  Hayiaond  L.  Biapllnghoff i  Profasaor  of  Aaronautlcal 
Dn^aaris^t  MITt 

*It  is  not  unreasonable  to  say  that  ve  have  reached  the  place  in 
the  development  of  a  nev  technique  vhere  it  is  no  lon^ar  necessary  to 
Justify  its  validity  or  usefulness.  It  is  true  that  thei'e  are  still 
unanswered  questions  and  inconsistencies,  but  the  advantages  to  be 
derived  from  the  application  of  power  spectral  analyses  to  the  gust 
problem  far  outweigh  these  short  comings.  The  principal  problem  to  be 
faced  now  is  how  the  process  can  be  placec.  more  efficiently  in  the 
hands  of  the  structural  designer.  This  is  indeed  the  most  difficult 
phase  of  the  problem.  In  its  present  state,  power  spectrum  analyses 
can  and  should  serve  as  on  adjunct  to  other  methods  proscribed  by 
regulatory  agencies,  and  one  vrould  expect  that  progressive  engineers, 
especially  those  dealing  with  unconventional  designs,  would  do  so. 

Dynamic  analyses  by  the  power  spectrum  api-roach  has  not  at  the  present 
writing  reached  the  state  v;here  it  can  be  applied  with  confidence  to 
the  art  of  airplane  design  as  a  completely  rational  and  independent 
process.  However,  its  present  value  to  engineers  in  comparing  the 
properties  of  one  aircraft  with  those  of  another  should  not  be  under- 
estimated.  It  has  unquestionably  reached  a  state  wilero  serious  consid¬ 
eration  should  be  given  to  incorporating  it  in  somo  form  in  government 
design  criteria. 

Three  places  are  suggested  where  effort  is  required  to  improve 
this  situation.  The  first  is  continued  effort  to  classify  the  atmosphere  - 
especially  under  conditions  of  severe  turbulence.  Secondly,  more  explicit 
data  on  the  range  of  turbulence  scales  to  be  encountered  is  needed.  Fin¬ 
ally  there  should  be  analyses  of  the  time  histories  of  older  aircraft  by 
means  of  which  the  power  spectral  dynamic  analyais,  bhou!^d  be  applied  to 
older  airplanes  whose  strength  and  service  life  are  well  known.  In 
this  way  reference  points  can  be  established  which  reflect  the  past 
experience  of  design  and  service  and  from  which  the  same  methods  of 
analysis  may  be  used  in  a  comparison  and  extrapolation  process  involving 
new  designs.  This  is  a  large  undertaking  but  perhaps  if  it  were  under¬ 
taken  by  the  builders  of  the  aircraft  in  question,  it  would  be  accennp- 
lished  in  reasonable  time". 


III.  Limitations  ot  Early  PSD  Analyses 


Perhaps  the  most  important  ooission  in  early  analyses  was  the 
pltohin^  degree  of  freedom  and  assooiated  short  period  damping}  eapeoi* 
cdly  for  unoonventional  configurations  involving  swept,  delta-type  or 
tail -loss  aiTangeaentB.  Also  the  effect  of  span -wise  variation  of 
gusts  la  important.  The  span-wise  variation  is  often  represented  by 


For  larger  values  of  b/L,  the  average  velocity 
progressively  reduced  particularly  at  the  high  fre- 


Quenoies  and  attenuation  in  wing  lifts  and  in  overall  airoroft  accelera 


tions  will  result. 


Not  always  so.  Actually,  the  effect  of  taking  the  span  wise  varia¬ 
tions  of  gust  intensity  into  account  is  to  reduce  the  mean  squared 
bending  looment  when  most  of  the  mass  is  in  the  fuselage  and  to  increase 
it  when  most  of  the  mass  is  in  the  wing. 

The  magnitude  of  the  effect  for  a  flexible  wing  appears  to  be 
greater  than  for  a  rigid  wing  where  the  analysis  is  based  on  the  first 
syama trial  free-free  bending  modes. 


IV.  Extensions  of  PSD  Analyses  and  Comparison  with  Flight  Results 

PSD  Anp-lyaifl  of  the  gost-load  problem  was  extended  by  Dloderlch 
(Ref.  7)  ^0  include  the  effect  of  lateral  veiriatioa  of  the  instanton- 
eoua  gust  intensity  on  the  aerodynamic  forces.  The  forces  obtained  in 
this  manner  vere  used  in  dynamic  analysis  of  rigid  and  flexible  air¬ 
planes  free  to  move  vertically,  in  pitch,  and  in  roll.  The  effect  of 
the  interaction  of  longitudinal,  normal,  and  lateral  gusts  on  the  wing 
stresses  were  also  considered.  .  ■ 

The  mean  square  values,  correlation  functions,  and  power  spectra 
of  soiaa  of  the  aerodynamic  forces  required  in  this  type  of  analysis 
were  calculated  for  one  special  correlation  function  of  the  atmospheric 
turbulence.  It  was  shown  that  if  the  span  is  relatively  large  compared 
with  the  integral  scale  of  turbulence,  the  mean  square  lift  and  root 
bending  nment  directly  duo  to  the  goat  are  substantially  reduced  when 
the  differences  in  instantaneous  intensity  of  the  turbulence  along  thq 
span  are  taken  into  account.  However,  if  the  motions  of  the  aircraft 
are  taken  into  account  the  mean  square  root  bending  moment  may  be  in¬ 
creased  as  a  result  of  these  differences. 

Also,  the  mean  square  pitching  moment  was  shown  to  be  substantially 
Increased  if  the  tail  length  is  relatively  large  compared  with  the  scale 
of  Uirbulence.  Finally,  the  wing  stress  due  to  longitudinal,  normal 
and  lateral  gusts  were  shown  to  be  statistically  independent  under 
certain  conditions.  (Ref  Fig.  1,  2,  3) 

The  effects  of  wing  flexibility  on  the  wing  strains  that  result 
from  gust  encountered  were  examined  by  Shufflebarger  (Fig.  4)  on  a 
power  spectral  basis.  Both  flight  results  and  analytical  results  were 
considered  for  a  four  engine  bomber  aircraft,  and  correlations  between 
the  measured  and  analytical  results  were  made.  The  effects  of  ■  wing 
flexibility  on  the  wing  strains  wera  moaourod  in  terms  of  amplification 
factors  based  on  the  ratio  of  the  strain  for  the  flexible  cohdition  to 
the  strain  for  the  "rigid"  condition,  and  results  wore  obtained  for 
four  spanwise  stations  of  the  wing.  The  effect  of  including  the  second 
symmetric  bending  node  in  addition  to  the  fundamental  bending  mode  in 
the  calculations  was  also  evaluated.  The  study  led  to  the  following 
results. 

1.  Comparison  of  the  measured  strain  spectra  for  the  flexible 
and  rigid  conditions  showed  strain  amplificationo  in  the  frequency  range 
below  0.5  cps,  especially  at  the  outboard  wing  stations,  that  ore 
attributed  to  the  effect  of  factors  such  as  spanwise  gust  variations  and 
aircraft  roHing  motions  rather  than  to  flexibility  effects. 

2,  The  measured  strain  amplification  factors  due  to  wing 
flexibility  at  a  station  near  the  wing  root  were  I.09  on  the  basis  of 
the  rms  strains  derived  fresn  distributions  of  strain  peaks.  The  ampli¬ 
fication  factors  decreased  with  each  successive  outboard  station  and 
then  increased  slightly  at  the  most  outboard  measuring  station. 


3.  1116  calculated  output  spectra  were  in  reasonable  agreeoeot 
with  tb«  oeasurdd  apeotra  except  at  lov  frequencies ,  where  the  diffoi*- 
enoea  are  attributed  to  oxtraneoui  faotoro.  At  the  outboard  atationa 

it  vaa  necessary  to  include  the  effect  of  the  second  syErnetrio  bending 
node  in  the  calculations  in  order  to  adequately  ooticate  the  effect  of 
flexibility  on  the  distribution  of  strain  peaks.  i 

•  ^  I 

4.  Both  calculations  and  neasurenents  indicated  that  the 
second  syjnnstric  bending  code  had  a  very  small  effect  on  the  nas  strain* 
but  that  for  strain  amplification  factors  detorminod  from  distributions 
of  strain  peaks,  the  second  symmetric  bending  mode  bad  a  pronounced 
effect  at  the  outboard  wing  stations. 

5.  Although  uncertainties  in  the  calculations  may  arise  ftraa 

a  number  of  factors;  these  results  Indicated  that  calculations  based  on 
PSD  nay  provide  reliable  estkates  of  the  magnitude  and  character  of 
wing  flexibility  effects  on  wing  strains  in  gusto  for  unsvept-ving  air* 
craft.  Also,  a  flight  investigation  was  made  by  Rhi'ne  on  a  large 
swept-wing  bomber  air^^lane  in  rougli  air  at  35*000  5*000  feet  to 

determine  the  effects  of  wing  flexibility  on  wing  bending  and  shear 
strains.  (Ref.  Fig.  5)  lu  order  to  eyaluate  the  overall  magnitude 

of  the  aeroelastlc  effects  on  the  strains  and  their  variation  with 
spanwise  location,  amplification  factors  defining  the  ratio  of  the 
strains  in  rough  air  to  the  strains  expected  for  a  'rigid"  and  'quasA*- 
rigid*  aircraft  were  determined.  The  results  obtained  indicated  that 
aeroelastlc  effects  are  rather  large,  jiurticularly  at  the  out-board 
stations.  The  effects  of  dyioamic  aeroelastlcity  Qp^.ear  to  increase 
the  strains  from  0  to  I70  percent  depending  upon  the  spanwise  station. 

On  the  other  hand,  the  relieving  effects  of  static  aeroelastlcity  appear 
to  reduce  the  strain  amplification  in  rough  air  by  a  significant  amount* 
The  response  of  swopt-wing  aircraft  in  rough  air  involves  a  number  of 
complications  not  present  in  the  case  of  unswept -v/ing  aircraft  due- 
princlpally  to  the  increased  importance  of  torsion  for  unswept-ving 
aircraft;  the  torsion  in  turn  results  in  significant  effects  on  both 
the  aircraft  dynamics  and  stability;  in  addition,  the  aircraft  vibratory 
modes  may  no  longer  be  approximated  by  simple  beam  bending  theory  but 
may  require  consideration  of  coupled  bending-torsion  modes.  Few  eiporl** 
mental  data  exist  on  the  character  and  magnitude  of  these  problems. 

Results  at  33,000  feet  were  compared  with  previous  results  at 
5,000  feet. and  comparison  showed  wing  trains,  on  the  average,  about 
20^  larger  at  the  higher  altitude.  Representative  values  of  the 
amplification  factors  varied  about  I.3  at  the  root  equations  to  about 
2.5  at  the  midspan  stations. 

yhe  analysis  of  the  strain  responses  of  a  large  swept-wing  airplana 
in  rough  air  by  Coleman,  Press  and  I-Ieadov/s  by  meaiiS  of  experimentally 
determined  frequency  response  functions  indicated  that  the  wing-b(5ndlllg 
and  shear  strain  responses  at  the  various  stations  are  amplified  by 
rather  large  amounts  because  of  the  responses  of  the  structure* 


I 


The  mount  of  amplifications  in  bending  strains  was  about  10  to 
20  percent  at  the  root  stations  but  increased  to  values  in  excess  of 
lOQ^  in  scxiifi  cases  at  the  midspan  stations.  The  shear  strains  shoved 
a  similar  pattern  across  the  airplane  span  but  also 'indicated  larger 
variations  between  front  and  rear  span  stations.  The  large  variations 
in  strain  responses  across  the  airplane  span  indicated  that  strain 
distributions  in  gusts  are  very  different  under  rough  air  loading 
conditions  than  under  the  usual  maneuver  loadings  and  warrant  detailed 
and  separate  consideration  in  design.  In  general  the  predominant 
source  of  strain  amplification  was  associated  with  the  excitation  of 
the  fundamental  wing-bending  mode.  However  at  the  outboard  statioao 
and  particularly  in  the  case  of  the  shear  strains,  significant  conta?i- 
butions  to  the  strains  arise  from  the  higher  syniaetrical  and  anti- 
symmotricol  vibration  inodes.  Thus,  the  effects  of  these  higher  modoo 
on  the  strains  may  also  have  to  be  considered  in  stress  oalculotioaie 
depending  on  the  degree  of  accuracy  rec^uired,  A  detailed  analysis  of, 
the  reliability  of  frequency-responso  function  estimates  obtained  by 
random-process  techniques,  particularly  as  affected  by  extraneous 
noise,  vrao  given.  The  effects  of  such  noises  in  giving  rise  to 
atic  errors  or  distortions  and  random  sampling  errors  were  explored  a 
and  results  of  general  applicability  obtained.  These  resiilts  wars 
also  applied  to  the  present  test  data  in  order  to  establish  their 
reliability  and  to  derive  adjustments  for  the  distortions.  The  Impart 
tont  result  obtained  is  the  indication  that  with  appropriate  praoau- 
tions  flight  tests  in  rough  sir  of  few  minutes  duration  may  bo  ufied 
to  obtain  reliable  estimates  of  airplane  frequency -re spouse  funotlonfl# 

« 

The  determination  of  aircraft  frequency -response  functions  for 
responses  to  atmospheric  turbulence  from  measurement  in  continucua 
rough  air  involved  a  relatively  new  application  of  random-prooeas 
"techniques.  A  general  analysis  of  the  reliability  of  such  freqiuenoy 
response  estimates  was  presented  and  methods  of  estimatii^  the  dlatoo?* 
tion  and  sampling  errors  were  developed  and  applied. 


V,  Evaluation  of  riothode  of  Calculating  Bending  Moment  Rospopsg  of 

Is^ Affl  e 


'Sura  to  arise  as  a  practical  Quastlon  whan  revitviog  and  gooparing 

analytical  raaulta  is:  ■Which  results  are  bettor?* 

Several  approxiinate  procedures  for  calculating  the  bending  mozaentj 
response  of  flexible  aircraft  to  continuous  isotropic  turbulence  wera. 
presented  and  evaluated  by  Bennett  (Reference  2).  Those  apiaroxinatol 
procedures  were  applied  to  a  slxaplified  aircraft  which  consisted  of  a 
uniform  beam  with  a  concentrated  fuselage  moss  at  the  center.  The 
conclusiC'ns  drawn  from  this  study  were: 

1.  The  force -summation  raethod  based  on  one  natural  bending 
mode  gives  very  good  results  compared  with  the  exact  solution. 

.  2.  The  matrix  method  based  on  five  stations  across  the  aemi- 

span  gives  very  good  results  compared  with  the  exact  solution. 

3.  The  mode -displacement  method  based  on  a  natural  raodo 
yields  inaccurate  results;  however  this  relatively  simple  method  can 
be  useful  in  trend  studies  involving  variations  in  wing  flexibility. 

4.  The  force -summation  method,  based  on  an  approximate 
parabolic  mode  and  the  Rayleigh  frequency,  loses  little  accui'acy  in 
mean-square  results.  Hov;ever,  the  approximated  natural  frequency  causes 
a  shift  in  the  fundamental  mode  resonant  peak;  consequently,  higher 
statistical  moments  will  be  in  appreciable  error. 

X 

''  5*  force -summation  method  based  on  an- approximate  para¬ 

bolic  mode  loses  little  accuracy  if  the  natural  frequency  is  kno-wn.  The 
inclusion  of  spanvrf.Be  variations  of  turbulence  results  in  a  decreased 
response  if  npst  of  the  aircraft  mass  is  in  the  fuselage  and  in  an 
increased  response  if  most  of  the  mass  is  in  the  wing.  Also  the  exclu¬ 
sion  of  spanwise  variations  of  turbulence  in  trend  studies  such  as 
trend  of  amplification  factor  -with  mass  ratio,  may  lead  to  erroneous  con¬ 
clusions. 


TI.  g^ety  flM  aaffty  natQrt.tQrMTrrsa^. 


It  lo  gtnartlly  aooepttd  that  guata  art  tha  dcrvjlnating  oauaa  of 
fatigua  vith  raapaot  to  tha  visga  of  transport  airplioiat,  vhila  tha 
dooinant  fatigua  dozoago  of  vings  of  flgbtar  or  trainar  planaa  is  dua 
to  caneinrer  loads* 

Tha  purposa  of  tha  applloatlon  of  statlstloal  mathoda  to  prohlaos 
of  trwiximum  loads  ajud  of  load  spaotra  is  to  obtain  by  extrapolation 
astixnatas  of  tha  probability  of  encountering  the  extrezaoly  hi^  design 
load  factors  of  vary  low  expected  frequency  of  occurrence  as  veil  as 
to  obtain  estimates  of  the  rata  of  accumulation,  under  random  loading 
.associated  vith  certain  load  spectra  or  fatigue  or  creep-damage .  It 
.should  be  obvious  that,  because  of  the  necessity  of  such  extrapolation, 
ithe  selection  of  a  probability  function  involves  more  than  curve  fitting* 
Unless  the  selected  probability  function  is  germane  to  tha  problem  and 
adequately  represents  the  inherent  statistical  variability  of  the  phenona* 
non  vhich  results  from  certain  basic  assumptions  concaming  its  origin, 
extrapolation  tovards  the  eactremss  (tails  of  tha  function)  vill  result  t 
.in  erroneous  predictions  vithin  this  range  of  variation,  vhich  is  just 
the  relevant  design  range.  It  appears  that  tha  probability  distribution 
of  marlmum  values  obtained  from  successive  samples  of  records  can  be  used. 
for  an  effective  estimate  of  the  required  extremes  for  design*  The  sta¬ 
tistical  theory  of  the  distribution  of  extreme  values  indicates  that  for 
initial  distributions  of  the  most  common  type  \diich  approach  xero  expon¬ 
entially,  a  limiting  form  exists  for  the  distribution  of  the  maxiinum 
values  of  large  samples*  The  use  of  this  asymptotio  form  should  there¬ 
fore  increase  the  reliability  of  the  estimate  in  comparison  to  direct 
^extrapolaUon  from  the  observed  initial  distxdbution* 

\  ,  -  ■  ■ 

The  character  of  tha  probability  function  of  extreme  gust-inten- 

aities,  with  the  aid  of  vhich  the  ultimate  load  factor  can  be  estimated, 
is  thus  ^oompletaly  deteimined  by  the  extremal  nature  of  the  phencoenon* 

It  is  independent  of  the  nature  of  the  underlying  distribution  of  gusts* 

If  .the  underlying  distribution  of  gusts  is  required  it  can  either 
be  obtained  by  statistical  interpretation  of  actual  .gust  records  or  by 
a  theoretical  study  of  the  effects  on  the  airplane  of  atiaosphsrio  turbtt- 
lonce  and  of  buffeting,  and  of  the  response  of  the  airplane  to  these 
random  disturbances  by  means  of  generalized  haxmonio  analyais  end  power 
spectra  techniques. 

*  •  ,.i 

The  design  gust  spectrum  sho\xld  extend  from  the. limit  load  downward* 
However,  the  small  probability  of  actual  occurrence,  during  the  opera¬ 
tional  life,  of  loads  in  the  vicinity  of  the  limit  load  makes  this  load 
region  rather  insignificant  vith  respect  to  damage  accumulation  over  this 
life.  The  specific  shape  of  the  spectrum  is  therefore  less  significant 
in  the  limit  load  region  than  in  the  region  of  higher  frequencies;  thus, 
discrepancies  between  actual  load  distribution  and  the  assumed  form  of 
the  spectrum  can  bo  tolerated  in  the  vicinity  of  the  limit  load.  The 


oonolujion  appearf,  th«refor«i  to  bo  justified  that  the  initial  fre(iuena7 
function  of  gusts  can  bo  fairly  veil  roprosantod  by  a  simple  exponential 
probability  function  over  the  entire  range*  The  saioe  can  also  be  dc»M 
for  nuu^ouvor  loads  for  fighters. 


VII  Ileaffuronect  and  Aasosscent  of  Bapeatod  Loads 

TUrther  li(^t  vas  thrown  on  tha  principal  problam— placing;  PSD 
techniques  in  the  hands  of  the  structural  designer—by  I'oneXy  (Rftf*  6)* 
He  outlined  the  General  principles  vital  to  the  xaaaaureaent  and  asaeu- 
icont  or  repeated  loads  as  follows  t 

!•  Lxtrcrwil  loads  that  oocaproaise  the  load  history  of  an  airplane 
should  be  treated  according  to  load  souroe*  Indications  are  that  there 
is  no  \inique  strain  relation  applicable  to  all  points  in  the  struotiire 
or  all  load  sources. 

2.  Tha  relativb  importance  of  the  different  load  sources  will 
vary  according  to  cd.asioQ  and  usage  of  ths  airplane. (Ilg.  7)* 

3*  It  is  mandatory  to  subdivide  the  load  histories  further  sooord- 
ing  to  flight  condition  or  a  soxu'ca  providing  more  hcncoganeous  infoma- 
tion.  Loads  and  their  sources t 


Load  Source 


Oust  Loads  cloud  and  wind -shear  turbulsnoe,  other  meteoro¬ 

logical  factors 

flight  conditions  -  climb,  cruise «  and  descent. 
In  the  case  of  climb  and  descent  although  the 
gust  experience  should  be  the  seme,  the  differ¬ 
ence  in  flight  speed  has  s  siguifioant  influ¬ 
ence  on  the  load  distributions 


Groiuid  Loads  air-to-ground  cycle,  landing  impact,  taxi  loads 


4*  On  th#  assumption  that  the  load  histories  of  futxtre  aircraft 
are  of  prime  interest,  the  aim  of  much  of  today's  work  is  aimed  st 
ooUeoting  data  on  ths  distrubanos  and  operating  conditions  rather  than 
oolleoting  masses  of  load  statistics  on  yesterday's  designs. 

Tbs  final  principle  guiding  research  in  connection  with  repeated 
loads  is  to  sxibdiTids  ths  work  into  studies  of  the  distrubenoe,  operating 
oonditions,  and  airplane  behavior.  The  stuiy  of  airplane  behavior  seems 
to  be  the  moat  formidable  task  of  tha  three,  since  It  involves  extensive 
instrumentation  and  flight  tests  under  controlled  conditions.  Other 
pertinent  considerations  arei 


1.  Tar  gusts,  the  standard  expression  for  the  pcwstr  speotrum, 
doesn't  oorer  the  larger  vare  lengths  Tery  veil. 

2.  The  selection  of  a  single  value  of  L  assumes  that  the  spectral 
shape  for  atmospheric  turbulence  is  independent  of  source,  an  assumption 
requiring  additional  verification  particularly  for*  the  .turbxilenoe  of 
convective  clouds. 

3.  Profiles  for  taxi  loads  in  the  form  of  power  spectrumd  indicate 
a  rapid  decrease  in  power  as  the  wavelengths  decrease.  A  major  diffi* 
culty  arises  from  the  fact  that  shook  struts  appear  to  have  a  threshold 
below  which  they  do  not  operate. 

4.  Por  couzoeroial  operations,  the  short -duration  lov->altitude 
flights  characterizing  feeder-line  operations,  gust  loads  consume  about 

of  the  airplane  fatigue  life.  For  high -altitude  let  transports 
6oe(  of  fstigus  life  can  be  ascribed  to  gusts.  (Ref.  o) 


VIII  Son0  Exporiences  in  the  Application  of  PSD  to  r>ancuver  acd 
132LL  Loads 

Application  of  PSD  to  maneuver  loads  obtained  on  jet  fighters  was 

IsTCfitiGAtad  by  Haiscri 

It  was  determined  that  the  maneuvering  load -factor  tixae  histories 
appeared  to  be  described  by  a  truncated  normal  distribution,  (fig.  8) 

The  power  spectral  densities  obtained  were  relatively  level  at 
frequencies  below  0.03  cycle  per  second  and  varied  Inversely  with  approx¬ 
imately  the  cube  of  the  frequency  at  the  higher  frequencies.  In  general, 
the  frequency  content  was  very  lov;  above  0.2  cycle  per  second.  The  load- 
factor  peak  distributions  were  estiir-ated  fairly  veil  from  the  spectrum 
analysie.  In  addition,  peak  load  data  obtained  during  service  opera¬ 
tions  of  fighter-typo  airplanes  with  flight  tine  totaling  about  24,000 
hours  were  examined  and  appeared  to  agree  reasonably  well  with  the  type 
of  equations  obtaineu  from  peak-load  distributions. 

Also,  the  frequency  content  of  some  airplane  response  quantities 
obtained  from  a  number  of  operational  training  flights  of  a  fighter 
airplane  was  presented  in  another  study  by  Hamer.  Power  spectral  densi¬ 
ties  of  normal  and  transverse  load  factor  and  pitching  acceleration 
wore  shovra  for  several  types  of  missions  normally  performed  by  the 
aircraft.  The  frequency  content,  which  is  described  by  the  spectrum, 
provides  information  which  is  useful  in  the  design  of  recording  and 
computing  equipment  for  analysing  maneuver-load  dara. 

V/hen  normalized  by  dividing  the  mean-square  value,  the  results 
indicated  that  except  for  some  differences  at  thelhigher  frequencies 
duo  to  the  effect  of  rough  air,  the  frequency  content  of  each  of -the 
airplane  response  quantities  v;as  similar  during  the  different  types 
of  mission  investigated. 

TI.0  normal -load -factor  data  for  the  different  types  of  missions 
exhibited  some  of  the  characteristics  of  a  Gaussian  random  process; 
therefore,  power  spectral  methods  were  used  in  analyzing  the  maneuver¬ 
load  factor  data.  Horraal -load -factor  peak  distributions  were  estimated 
to  a  reasonable  degree  of  accuracy  from  this  spectrum  analysis.  Peak 
distributions  for  transverse  load  factor  and  pitching  acceleration  were 
not  determined  from  the  power  spectrum  because  in  most  typos  of  missions 
the  quantities  did  not  appeaf  to  have  the  charaoteristios  of  a  Gaussian- 
rondom  process. 


Morris  (Ref.  21)  conducted  a  study  of  the  response  of  a  light  air¬ 
craft  to  roughness  of  unpaved  nmirmys;  and  found  that  mean-v/quare  accel¬ 
eration  response  increased  approximately  linearly  up  to  about  2S  or  30 
mph.^,  but  increased  more  rapidly  than  the  linear  rate  at  higher  speeds. 
Transfer  functions  computed  from  the  air  craft  acceleration  response  to 
the  runv/ays  v/ere  in  relatively  good  agreement  at  20  mph.  As  regards  to 
the  peak  amplitude  and  frequency  at  which  it  occurred,  apparent  lack  of 
agreement  at  1^0  mph  may  have  resultea  from  nonlinearities  in  the  system 
and /or  for  small  a  statistical  sample  for  determining  a  transfer  function. 


IX  A5Ri»PS~l]L  Experiences  in  The  Application  of  Powora  Spectral  l*bthods 
to  Alrplano  Dasign  under  Prolect  116?  Structural  IhBirn  Criteria 

Stulv  of  World  Vide  Turbulence  Date  (noferonco  2?) 

Nev  York  University  is  oaking  an  analysis  of  B-66B  low  level  turbulence 
data  to  determine  spectral  variance  and  spectral  slope  behavior  with  respect 
to  four  physical  paremstorj  Mean  wind  speed,  atmospheric  .stability,  terrain 
and  height}  and  is  re-examining  the  basic  approach  to  the  gust  load  design 
problem.  Current  data  indicates  that  the  normalized  spectral  shapes  are 
sensibly  independent  of  wind  speed  and  stability.  Additional  analyses  are 
necessary  to  confirm  these  spectral  shape  results. 

Heating  Effects  on  Gust  Response  (Reference  26) 

The  response  of  a  high-speed  flexible  vehicle  to  atmospheric  turbulence 
was  investigated  for  flight  at  an  altitude  of  1000  ft.  The  Mach  number 
range  was  from  0.7  to  6.0.  The  response  of  the  structure  at  nominal  room 
temperatures  was  compared  with  the  response  of  the  structure  when  ae rodynamically 
heated  to  a  uniform  temperature.  (Reference  Figure  9)  The  response  to  discrete 
guets  was  compared  with  the  response  to  continuous  random  turbulence  (computed 
by  power  spectral  densities),  and  the  effect  of  flexibility  was  investigated. 

Both  ccthoda  used  oscillatory  pressure  distributions.  The  results  have  not 
been  tested  and  are  applicable  to  similar  configurations  only.  These  results 
indicate  that:  (a)  Rms  acceleration  and  bending  moment  responses  do  not 
necessarily  increase  with  speed  continuously,  (b)  Maximum  rms  ^nerally 
occur  at  between  Mach  4  and  ^iach  5,  except  for  fuselage  bending  moment 
which  is  slightly  higher  at  r-lach  6,  (c)  One  degree  of  freedom  underestimates 
the  rms  acceleration  and  bending  moment  responses  with  the  exception  of  the 
fuselage  bending  moment  in  the  Mach  number  range  between  Mach  3  Mach  5» 

(d)  The  number  of  degrees  of  freedom  have  a  considerable  effect  on  the  rms 
response.  An  extreme  example  is  the  rms  nose  acceleration  at  I^ch  5*  The 
five  degrees  of  freedom  rms  response  is  approximately  33^%  higher  than  pro- 
dieted  by  one  degree  of  freedom. 

V/hereas,  the  discrete  gust  method  has  no  means  of  predicting  a  change 
in  the  number  of  peaks  exceeding  the  zero  level  due  ■co  heating,  the  con¬ 
tinuous  gust  method  does  predict  a  change  in  the  number  of  peaks  exceeding 
the  zero  level  through  the  change  in  the  characteristic  frequency.  The 
rms  response  level  of  nose  acceleration  increases  by  approximately  62% 
duo  to  heating.  The  ratio  of  the  fuselage  bending  moment  stress  to  the 
ultimate  stress  increases  by  approximately  75%  due  to  heating.  Over  one- 
half  of  the  latter  increase  is  due  to  the  reduction  in  ultimate  stress 
produced  by  heating. 

Optimum  Fatigue  Spectra  (Reference  23) 

Two  types  of  power  spectral  density  analyses  of  taxi  loads  were  made. 

The  first  was  an  analog  study  in  which  the  complete  non-linear  system  of 
gear  equations  was  considered.  Runways  were  generated  and  the  response 
of  the  system  in  the  time  domain  to  taxi  over  these  runways  was  produced. 
Cumulative  occurrences  of  incremental  c.q.  acceleration  levels  were  then 
obtained  by  counting  the  number  of  crossings  with  positive  slope  of  the 


various  levels.  The  second  typo  of  analysis  coiqpri^'Qd  a  ^digital  study* 

It  was  the  aim  of  the  analysis  to  ascertain  whether  it  was  possible  to 
establish  the  ouzculativo  occurroncoe  of  incremental  c*q.  acceleration  by 
conventional  spectral  analysis  in  the  frequency  docBin  with  a  suitable 
llneorized  representation  on  the  main  and  nose  gear  forces.  The  airplane 
configuration  used  was  the  DC-7o  at  143 gross  weight.  A  single 
taxi  velocity  of  22  ft  per  second  was  considered  throughout.  The  results 
of  the  digital  analysis  wore  found  to  bo  powerfully  effeotod  by  the  line¬ 
arized  representation  of  the  oleo  dampings*  Linearization  of  the  oleo 
dajg>lng  requires  the  selection  of  an  amplitude  and  frequency.  In  AGAKD 
heport  119 •  it  is  suggested  that  the  basis  of  selection  should  be  a 
dominating  strut  amplitude  and  frequency.  A  new  procedure  was  arrived 
as  to  obtain  a  more  specific  criteria  for  linearization  of  the  oleo 
damping,  having  due  regard  to  the  strolce  response  at  all  frequencies. 

By  way  of  comparison  the  non-linoor  results  exhibit  a  definite  bias  toward 
the  negative,  or  upward  accelerations.  For  negative  c.g.  accelerations  the 
agreenunt  between  linear  and  non-linear  analysis  was  satisfactory.  The 
most  promising  procedure  seems  to  be  to  consider  the  non-linear  gear  with 
a  stationary  random  input.  The  non-linoar  gear  analyses  also  show  that 
the  resonances  associated  with  the  unsprung  gear  masses  can  have  an  appraol- 
able  effect  on  the  cumulative  occurrences  of  c.g.  accelerations  unless 
suitably  filtered.  These  resonances  ordinarily  lie  at  frequencies  beyond 
the  range  of  definition  of  the  runway  roughness  power  spectra  of  NACA  TN 
4303*  Oscillograph  records  obtained  during  airplane  taxi  conditions  have 
shown  high  frequency  content,  indicating  excitation  of  mode  lying  outside 
the  range  of  definition.  It  would  therefore  seem  desirable  to  extend  the 
bandwidth  of  the  spectra  of  NACA  TN  43^3  'to  higher  frequencies*  It  should 
be  re-emphasized  that  the  above  results  were  based  on  one  taxi  velocity. 
There  is  c  lack  of  experimental  data  on  the  distribution  of  taxi  Velocity. 
Moreover,  the  direct  application  of  runway  data  in  spectral  form  to  fatigue 
studies  is  complicated  by  the  non-linear  transmission  characteristics  of 
landing  gears  in  current  usage.  Other  phases  of  taxi  loads  such  as  tumingi 
pivoting,  braking  and  towing  should  be  considered  in  fatigue  analysis. 
Unfortunately,  the  only  data  available  on  these  afe  from  specific  flight 
test  programs  conducted  to  verify  design  criteria,  in  the  traditional 
sensei  not  fatigue  criteria*  Much  analytical  insight  has  been  gained 
by  investigations  to  date,  however. 

B-66  Low-Level  Gust  Study  (Reference  3I) 

Statistical  design,  data  gathering,  instrumentation,  data  reduction, 
ooxcputation,  and  presentation  were  accomplished.  Technical  analysis  of 
the  data  was  performed,  i.e.,  determination  of  the  scale  of  turbulonot, 
numerical  filtering,  integration,  Tesolution  and  spectral  estimation. 
Comparison  of  current  Mil  Specs  on  low-altitude  gusts  and  response 
equations  with  sampled  data  results,  recommendations  for  future  analysis 
and  flight  programs.  Comparison  with  Mil  Spec  8866  euad  NACA  TN  4332 
showed  that  sac^ile  distributions  of  standard  deviation  of  low-level  gusts 
differ  significantly,  equations  in  specifications  overestimate  the  number 
of  roost  frequently  occurring  gust  peaks,  which  are  of  interest  in  repeated 
loads,  and  underestimate  the  number  of  large  gust  peaks.  The  technioal 
analysis  further  showed  thatJ  Extermal  procedures  by  Oumbal  are  very 


useful  in  relating  expected  gust  peaks  to  miles  of  flight.  Isotropic 
turbulence  approximations  are  excellent  fits  to  auto^correlation  and 
power  spectrum  of  low-level  gusts.  The  Chi-square  distribution  is  a 
significant  fit  to  sample  distributions.  Sample  derived  equivalent  gust 
velocities  are  good  (not  tested)  estimates  of  observed  gust  velocities. 
Cross-spectral  power  between  sample  gust  velocity  components  in  negligible 
phase  is  uncorrelated]  amplitude  is  uncorrelated.  Transfer  functions  by 
cross  spectral  procedures  are  good  (not  tested)  estimates  of  analytical 
transfer  functions.  Statistical  techniques  such  as,  cross  correlation, 
sampling,  regression,  product  moment  correlation  coefficients,  etc.] 
are  useful  tools  in  separating  gust  and  maneuver  responses;  also  the 
distribution  of  gust  lengths  is  useful  in  separating*  gust  and  maneuver 
responses.  Plans  are  currently  being  made  to  gather  similar  data  for 
the  susziar  months  to  arrive  at  a  more  comprehensive  gust  model. 


X  RolatlnR  Peak  Cotmt  Data  to  RJ’S  Data 

In  ragarda  to  techniquea  for  relating  peak  count  data  to  naa  data 
for  lov-l07«l  flight,  Noteai  (Reference  24)  concludes  that  these  two 
types  of  data  (peak  count  and  rms)  can  be  related  in  a  way  which 
enables  a  siniple  application  to  the  specific  problem  at  hand.  Data 
from  'Notesa,  Analyses  of  Turbulence  Data  I'basured  in  Flight  at 
Altitudes  up  to  1600  ft  above  Three  Different  Types  of  Terrain,* 
and  CAL  Report  Nr  TE-1215-F-1,  Foburary  1939  WADD  TR  6O-303, 

*£-66B  Low  Level  Gust  Study,'  was  analyzed  and  indicates  that  on 
the  average  a  gust  of  extreme  mgnitude  is  usually  associated  with 
a  nearly  homogeneous  patch  of  turbulence  having  an  rma  of  one-fourth 
the  extreme  value  in  an  extent  of  ten  miles.  Thus,  it  is  reasonable 
to  assume  that  in  every  flight  there  will  be  at  least  one  patch  of 
turbulence  (10  miles)  which  is  associated  with  the  peak  gust  encountered 
and  will  have  an  rma  gust  velocity  equal  to  one-fourth  of  this  peak 
value.  One  peak  exceeds  four  rT every  ten  miles.  This  is  based  on  40 
data  runs  (from  B-66  data)  wherein  the  peak  variation  for  35 
runs  was  between  3*4  ^  to  3*0  ^  • 

Considerable  work  remains  to  be  done  so  that  a  workable  model  can 
be  obtained  which  will  enable  one  to  take  into  account  the  relationship 
of  the  intensity  levels  among  adjacent  turbulence  patches. 

Existing  turbulence  models  such  as  NACA  TR  1272  do  not  provide  data 
on  these  interrelationships  since  peak  count  data  for  many  flights  is 
lumiped  together  into  one  probability  distribution.  (Figure  9) 

This  can  be  done  for  design  purposes  to  some  extent  if  one  applies 
the  probabilities  to  missions  rather  than  to  individual  ten  mile  patches. 
In  other  words,  one  should  say  one  out  of  seven  missions  will  be  performed 
In  turbulence  of  five  CT  rather  than  one  out  of  seven  patches  in  a  given 
mission. 


Rooognlzing  that  the  reluotance  to  uae  power  apectral  oothoda  la 
laoatly  due  to  the  complezlt/  of  the  analytical  techniquaa»  and  the 
fact  that  the  techniquea  of  atatiatical  analyaia  are  often  unfanillar 
to  deal^ra.  It  la  Inoumbent  upon  the  Air  Force  to  puraue  inyaatlgatlon 
aicud  toward  simplifying  and  applying  power  apectral  techniquea*  Toward 
thla  end,  Applied  Research  Ohjectlve  Nr  35A53  vaa  eatabllahed. 


Loaui.  gag  roBDoncP  rower  jppcvrua  oi  hp  AaTmeoa 

Vehicles 

To  determine  the  spectral  shape  and  Intensity  of  the  Input  and 
response  power  spectrum  of  advanced  vehicle  structures  result¬ 
ing  from  gusts,  ground  winds,  wind  shears,  pilot  Inputs,  auto¬ 
matic  guidance,  and  maneuvers.  In  order  to  apply  these  spectra 
and  root-maan  square  distributions  In  structural  load  calcula¬ 
tions  In  which  the  effect  of  vehicle  motions  and  flexibility 
are  Included  In  the  determination  of  vehicle  transfer  functions 
which  can  be  converted  to  real  time  differential  equations 
solvable  on  high-speed  coc^uters. 

PRESENT  STATUS!  The  spectral  shape  of  gusts  has  been  adequately  determined. 

Work  Is  continuing  on  the  variation  of  Intensity  with  geographic 
and  meteorological  factors  which  Influence  gusts*  Input  power 
spectrum  £iid  response  power  spectrum  are  being  measured  for 
groiind  winds,  wind  shears,  pilot  Inputs,  and  maneuvers.  No 
methods  are  currently  available  for  converting  available 
data  from  Its  present  form  Into  a  form  appropriate  for  spectral- 
type  calculations. 

TECHNICAL  POSSIBILITIES »  Further  analytical  Investigations  to  determine 
methods  for  re-evaluatlng  currently  available  data  on  winds, 
wind  shears,  pilot  Inputs,  automatic  guidance  systems,  and 
maneuvers  from  their  present  form  Into  a  form  appropriate 
for  spectral  calculations  as  has  been  done  for  gusts  In 
NACA  TR  1272.  Further  analytical  Investigations  to  develop 
a  mathematical  model  representative  of  the  spectral  shape  of 
wind  shear,  pilot,  maneuver,  automatic  guidance  system  Input. 
Further  analytical  Investigations  to  determine  goodness  of 
fit  criterion  for  sac^le  distributions,  of  this  data.  Further 
analytical  Investigations  to  determine  extremal  procedures  In 
relating  expected  peaks  to  miles  of  flight  for  Inputs  other 
than  gusts.  Further  Investigation  to  determine  the  significance 
of  cross-spectral  power  between  sample  wind  and  wind  shear  of 
cross-spectral  power  between  sajqple  wind  and  wind  shear  ooa>- 
ponents.  Further  investigation  to  determine  how  well  transfer 
functions*  Further  Investigation  to  determine  how  useful 
statistical  technique  such  as  cross-correlation,  sampling 
regression,  product  moment  correlation,  coefficients ,  etc*, 
are  In  separating  responses  from  the  aforementioned  inputs* 


I 
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Neither  ttiia  rapcrt  nor  th«  appliad  raaaarob  obJaotlTa  praauaoa  to 
anooopaaa  all  mv  eoneapta  vhiob  ii  auooaaafully  porauad,  vould  aahaAoa 
our  knovladga  and  adranca  tha-atata-of-tba*art. 

In  faoti  it  la  hopad  that  thla  raport  aanraa  io  diaoorar  and  atlsulata 
inraatigation  and  praaanta  aoaa  aithoda  of  daaonatration  vhioh  tand  to  laad 
tha  raadar  to  InTaatlgata  furthar. 
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APiliNUDC  II 


INrHOlXJCTlQN  TO  iiPhiCTJiAL  AUkUblH 


Sp«ctr.l  Dooalty  la  nothimi  iww. 
confrontad  with  tba  uaual  onglneoring  ^oblaa  of 
trying  to  diacovwr  a  pattarn  which  couplad  with  ingenuity  will  anah)* 

m.t  obyJom  form  of  ropotltlon  la  numloa  around  a 
of  «»pra»«  all  pradlctlona  In  tarma  of  tha  framuork 

P®^todlclty  or  irequency  and  aaplltude, 

y  =  a  sin  Wt 


analyj“a!““‘'  °t  ‘>.r»nlo 

„„h  ««tual  raoorda  of  a  contlnuoua  dl.turb.noa 

of  «  ««  dealing  with  tha  problem 

analyalo  of  a  random  tlnu  bletory 
Which  la  neither  periodic  (repetitive)  nor  decaying  with  tlna.  ^ 

do.lrJ'‘for^““  ‘.‘J''  =>>0100  but  to  OYorconu  tbe  natural 

I  •omotblng  nailed  down,”  and  to  apeak  tbe  language  of  pro- 
babillty  Instead  of  the  language  of  certainty  a 

Jr  randora,  It  rauat  have  elem^nta  of  uncertainty 

and  all  one  can  hope  to  determine  la  the  degree  of  uncertainty.  ^ 

.anted?  ““  "tbtlatlcal  terma  will  now  be  pro- 


STATlirriCAL  Tliiii-OWOiXXJY 


htAN  OH  AVmUKi 

Lut  t  Im)  a  raal  variable,  and  y  <•  f(T)  be  any  function  given 
on  an  interval  (a,b).  Then  7*,  the  average  of  y  with  respect  to  t  for 
the  interval  a,b  is  defined  by  ^ 

^  J  ^  j  t  Lt  }  dt 

^  ^  a 

(  b-a)  ^  ^ 


or 


a 


The  gaomatric  interpretation  nukes  it  easy  to  remember  the  definitioni 


for  positive  y(t) 


for  negative  y(t) 


bTATlsriCAL  Ttifu-lIhOUXiY 


hOOT-i-tAl^-iiv*UAWi  (hMi) 


Wo  (loflno  y,  tbo  root  totiaQ  square  (rma)  value  of  y  with  respect  to  t 
for  the  Interval  Stb  by  the  equation* 


Thus  the  iiiu  y  is  the  square  root  of  the  averut^e  of  the  square  of  y,  and 
is  so  defined  that 


This  definition,  also  beco*a3s  cluai'er  with  an  example  t 


Consider 


on  the  interval 


0,1.  Then 

100 

im+f 


so  that  in  this  case 

And  if  We  taRe  special  values  of  m. 


FOR  ^ 

4 

1 

0.ZS 

lEH 

0.224 

3.333 

5.773 

MU 

• 

o 

o 

The  extreme  values  of  y(t) 

=  |Ot*^on 

0,1  are  0  u 

iid  10.  This 

illustrates 

that  the  rms  value  of  y  always  lies  between  the  extreme  values  of  y(tj 
but  may  lie  quite  close  to  either. 


irrATIbl'ICAL  ThJu-U^OLCX;/ 


Ai/roooiuuUTioh  j-ui^riui^ 


Thu  uattleut  quantitativu  way  to  judge  the  corralatlon  or  fit  batvaao 
two  things  la  by  thalr  product.  In  guitural,  iurgu  products  man  good 
corrulatlon  aiid  small  products  moan  poor  corialatlon.  Whan  tba  'things* 
aiM  dependant  on  tim ,  say  f^tj  and  fit).  could  ba  aqual  to  fj,  t^ 'f*  T* 

ono  takas  tha  avaraga  px'oduct  ovur  all  tha  time  Involved*  Thus  tha 
autocorrelation  of  a  tlnu  history  Is  <m 


tLbM.WIl>  Ct'  ThAt»:>KQlu-ATION  THLOtfif 


lutlurniDti  to  tho  probluu  cf  tho  aoipll tuda-f ruquancy  aaalyaia  of  a 
tlDu-hlatory ;  it  may  Ui  elthor  periodic  or  not.  If  atrictly  periodic 
it  CdQ  usually  be  oxproasad  aa  a  Fourier  ::>uriea. 

% 

(Notice  how  we  now  have  related  tiur^,  t,  to  frequency,  CO  ,  and  vice  veraa.) 

If  f(  t)  la  not  atrictly  periodic,  it  luiy  either  aettle  to  aouu  constant 
value  as  tiiiu  on  or  it  nay  not. 

If  f(t)  ia  non-periodic  but  decays  to  zero  as  time  incruaaaa  it  can 
(usually)  be  expi'esaod  by  a  alit^it  generalization  of  Fourier  iieriea  using 
continuous  rather  than  discrete  iraqueuciea. 


LLfc:]-fcJ/Ib  OF  TKAKSFOhJ-ATION  TlUiOilY 


Tbla  traixafui-uiitioo  from  a  tluu  fuQctioo  to  a  froijuancy  functloo  ia 
ezpraaaad  by  tha  Fourier  intogral 

G(jw)=f"F(t)ej“^clt 

-CO 

For  any  function  F(t)*0,t  >0,  the  limitB  of  integration  can  be 
changed  and  tho  t  rune  form  i-ustated  oa  , 

0(jw)= 

maicing  use  of  the  Euler  rulutionshipi 

J(4t  . 


ueparuted  into  components 


M  f  J^I^Jcan  bu  ueparaiea  inio  components 

A w  S  Q  [G  (i  w)]  *  J  cos  u)t  clt 


0 
oa 


B«s  l[GCiw)]>f  J  F(t)  sin  wt  dt 
Q(jca)» 


In  polur  coordinates 


whe  re 


*1^ 

U) 


0  s  OJtC  tO/A. 


(•) 


The  practical  problem  ie  to  determine  a  mrithemtical  representation  for 
such  that  the  integrals. ( A)  (5)  evaluated.  Again,  the  question 

as  to  whether  J (t) 6 ^ exists  can  be  answered  aff irimtlvely  for 
any  F(t)  describing  a  physically  realizable  motion^ 

More  pertinent  is  tho  queatlom  Can  the  function  be  approximated  by 
F,(t)  such  that 


U) 

(v 

(^) 
(5) 

((J 

(7) 

(8) 
6(jca) 


0) 


some 


ELt.^tWIb  OF  ThA^SFOH^UTICW  THhOHY 


nearly  «o,  a/tor  aoa 

flnlta  tlnu  T.  If  this  coaditlon  cannot  be  nat  It  la  atlll  Doaalble 
to  app^unate  (2)  be  breaking  up  F(  t)  Into  parta  . 

're  fp(t)  ia  fluch  that  the  aecond  Integral  on  tho  right  can  be 


where  F^2^■t)  ^a  auch  that 
formally  evaluated 


Tbd  udXl  logical  quaatloo  iai  *If  thaau  traottforms  ar«  ao  easy 
to  avaluate,  why  did  wa  akip  ovur  that  partT* 

So  to  that  aod,  ao  iotuitiva  dascriptioo  which  is  consiatoot  with 
tha  laval  of  thia  raport  ia thown  below.  Thia  approach  waa  aalected  to 
avoid  the  'foraat-for-tha-traea*  eoigma. 


COi-lPl/rATICfiAL  TLCHNI^UC^  KOH  KOlftHIhlH  THANiiKOhi-V.TIC^J 


Given  a  function  (of  tluu).  fCt),  Its  Fourier  Transforo  F((a))  (•  function 
of  frutjuency)  is  defined  • 

wbsrs  ■  -1  and  e  Is  the  exf^onential  functions  We  shall  use  the  identity 

e*J^-COS  cat  -j  sm  cat 

and  shall  consider  sine  and  Cosine  Fourier  Transforms 

4«s  «*)«»'  at 

For  much  of  our  dlscusslout  we  consider  only  the  Cosine  Transforms  since 
a  parallel  discussion  would  apply  to  the  bine  Transform.  Furthermore, 
since  mutheoutlcally  the  tluu  and  frequency  domains  are  Interohantjeable , 
we  shall  sometimes  refer  to  F(O0)  as  the  bpectrum  In  order  to  emphasize 
that  It  Is  a  function  of  frequency.  This  usage  arose  because  the  motivation 
for  this  study  came  from  the  problem  of  obtaining  PbD  from  Autocorrelation 
Functions . 

Let  us  assuiifi  that  f(t)  Is  an  even  function)  l.e.,  f(  t)  ■  f(>t).  We  consider 
each  frequency,  and  aslc  how  much  relative  contribution  a  cosine  of  ttiat 
frequency  mokes.  We  consider  f(t)  to  be  made  up  of  a  number  of  component 
frequencies  (actually  a  contlnum  of  frequencies),  and  we  wish  to  determine 
the  Illative  amptltudes  associated  with  the  frequencies. 

The  way  In  which  the  cnlne  trunsform  answers  this  question  may  be  described 
us  shown  In  Figure  1.  Given  f(t),  we  choose  soioj  f requancy ,CJ ,  and  compare 
CO^  cot,  with  f(t).  In  order  to  make  the  comparison,  we  form  the  product  f(t) 

COi  cot  .  This  product  will  be  positive  when  f(t)  and  C05  Cjut  are  both 
positive  or  both  negative;  and  Its  values,  positive  or  negative,  give  an 
instuntnneoui  picture  of  whether  f(t)  andCOS  wtagixie.  But  we  want  to  compare 
f(t)  andCOi  Cot  throughout  their  whole  range,  ao  wo  integrate  the  product 
f(t)C0S  cot,  .  As  CO  is  varied,  this  process,  which  gives  a  step-by-step 
interpretation  of  the  definition  of  the  Cosine  Transform,  will  give  a 
relative  measure  of  how  much  contribution  each  frequency  makes  to  f(t). 

When  we  consider  the  process  otep-by-atep  in  this  fashion,  we  shall  refer 


Koi'ui  chu  iii:i idiituitcoua  ii«juiiur«)  of  cos  COi'tl 


InUif^'utb  lu  obtulu  the  lotui  iiuabul^  of  u^rbeunli I-  arca  *A 


A « {  (t)  cos  to^  dt  s  f  (la,) 

jt: 


This  is  thb  vdluo  of  Ihu  yourlsr  ICoslno)  truiisform  ul 
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binco  Hi)  auJCCfttjtare  oven,  l!\-(jO|)  =  F((jJv}. 
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Now  to  tie  tut{<*thor  aoim  of  the  topics  just  discuaaudt 

Thun  far,  in  our  oflort  to  uotelop  an  Hwpl I tuda- frequency  analyala 
for  a  random  tinw  history,  wo  rmvei 

1)  Dlaouaaed  the  autocorrelation  function  which  describea  the 
avera^  value  of  the  product  between  valuua  ol  the  record  a  specified 
tinu  interval  apart  aa  a  function  of  this  tim  interval  itaelfi  and 
therefore  givea  the  desired  amplitude  characterization  of  the  random- 
tinu  history ,  f( t) . 

2)  We  have  diacusaed  the  Fourier  Integral  Transformation  which 
provides  a  frequency  characterization  of  the  randoin- 1 iiaa  history. 

Wo  now  should  be  able  to  analyse  input- transfer- response  for 
random  functlonu  by  ejiprusaing  input  and  roaponae  as  Fourier  Integrals 
of  the  autocorrulutiona  of  tiieir  respective  tinu  histories.  In  any 
case,  we  can  impleirunt  this  "tongue-twister*  with  the  aid  of  the  follow¬ 
ing  outline  of  the  actual  nuchanics  normally  involved i 

THE  fiiCHANICii  OF  COllPUriNO  THE  POWEH  SHiCTHA 

1.  Hecording  time  sampled  input  data  on  magnetic  tape. 

2.  Feeding  this  dota  into  a  large  high-speed  computer. 

3.  Computing  the  input  spectra  from  this  data. 

4.  Computing  the  transfer  function  for  the  system. 

(If  the  input  and  output  spectra  are  known  the  transfer 
function  is 


3.  Hecording  the  output. 

6.  Computing  the  response  spectra. 


Tit.  TUtCbTf  PtiOCi^DUltL  FOh  COl-lWt^  PO»^li  bPtCTllA, 


For  a  f  ini  to  .^cord  lon^th  aumpled  at  t,  tho  following  atapa  may 
ba  uBod  in  computing  gust  spectra i  (although  other  nuthods  are  available) 


1.  Calculate  the  ouan,  y(t), 


whura  N  -  ^of  values  of  the  random 
variable  taken  at  equal 
units  of  tinui^t,  apart. 


2.  Calculate  the  autocorrelation  coefficients, 

where  M=#of  lags  or  uniformly 

spaced  points  over  the  frequency 
range  at  which  estimates  are 
desired 


-*  V  A  WM  w  ww  *  A  *  V  jav  M  was  I 

I  Vt^UlPS 


3*  Calculate  the  initial  eatiiaitea  of  power  spectra 


l-A- 

""  Calculato  • 


k*  Calculate  final  oatimates  of  power  spectra 


Je  *  0.5A  Lt  ♦  0.4 1  L, 

0.23  6^.,  +0.54  0.23  L 

O.qioL^^^i  +0.54  Lyri 

Calculate^^4*/'^Ql^®^  of 


>lt+f 


Thoao  power  apoctral  calculations  my  ba  quickly  obtained  by 
high-apeed  digital  computert 


POWLH  0>'  A  JtAWDO.-l  DlbrU]\bANU. 


Tho  pur;^otitf  of  thlu  ouctlon  la  to  llluatrato  tte  practical  aapacta 
of  tliu  Ol'uiVilajtitiOUOll  tupicu. 

If  a  ruiiduio  tiiiu  hlatory  la  conaldcrod  aa  boin^  analotjoua  to  current 
flowluti  ttiruU{^i  a  unit  ruaiatunce,  an  important  charac teriatic  la  the  power 
or  anurt^  conluinod  in  f(t).  Thia  powor  may  te  conceived  aa  composed  of 
cootributiona  from  various  frequencies.  One  speaks  of  the  density  of  power 
or  energy  as  the  limit  of  the  power  in  a  frequency  band  as  the  bandwidth 
approaches  z.ero.  This  plot  of  power  against  frequency  is  appropriately 
called  the  PbD  curve  or  power  spectrum,  by  knowledge  of  tho  transfer 
system,  the  input  power  spectrum  can  usually  be  converted  to  the  output 
or  response  power  spectruia^  aa  shown  above. 

In  order  to  develop  a  frequency  representation  which  would  be  applicable 
to  continuing  disturbances,  the  theory  of  random  processes  makes  use  of  the 
concept  of  a  atationury  rundoia  process.  The  basic  assumption  charaoterizin/ 
a  Btationury  random  procesu  is  that  the  underlying  mechanism  which  gives 
rise  to  the  disturbance  are  invariant  with  time  and  statistical  prediction 
becoioes  possible. 

Kor  tho  cuBo  of  a  stationary  random  function  of  time  ytt),  the 
moan  suuaivi .  is  defined  by  tho  following  equationi 

The  lamn  squai'e  will  usually  exist  and  represent  a  measure  of  disturbance 
intensity,  binco  y^(t)  is  a  quadratic  function  of  y(t),  it  has  been 
termed  the  "average  power"  of  y^Ct)  in  analogy  to  electrical  power  which 
is  proportional  to  the  square  of  the  current.  The  function  ylt)  is  con¬ 
sidered  to  bo  composed  of  a  infinite  number  of  einueoidal  components  with 
circular  frequencies ,  Co  •  between  0  and  infinity.  The  portion  of  y^tj 
arising  from  components  having  frequencies  Potween  CO  andCO  +  clwlB  denoted 
The  function  )  is  called  tho  BbD  function,  and  has  the 

property  that 

f  (aO  do) 

The  PSD  function  is  therefore!  ■  ^  I  ^ 

Equation  (3)  nay  be  used  to  evaluate  the  PSD  function  from  observed  data. 
However,  in  practice,  the  I’SD  function  may  be  more  conveniently  determined 
using  a  related  function,  the  Autocorrelation  function,  H(  ),  defined  by 


Thu  uuLucoiTdlutluii  i'liuctiuu  hud  thu  tr Icul  property 

uDd  lb  ruclprucully  iciutud  to  thu  I'uactluu  by  thu  Kouritfr  Coalao 
TruiibTui'u^tloit  lit  iho  iolluwlii(/  Uiuuiiuri 


{ 


RC'f)'irSC<4co»wy<*<o 

=  Vir  h’^R('d«s  ci'f  d.'f 

.a*  *.  »  »•..  »  A> t.  J ^  1  .  . 


Thlu  duflnltlon  ol'  thu  ioU  function  ib  conulbtuUt  with  tho  piucodlUfS 


defiiiltiou  of  u4Uution 
for  Kourlur  Trunaf ui'uuitlun. 


t>cu  Api>undlx  -  Computbtlonul  fuchni^uus 


When  a  linear  systum  la  uxpubuu  to  an  Input  varying  in  a  raudota  laannar 
with  tluiu,  thu  probability  dlatributiun  uf  thu  ayatum  output  y  can 
frequently  be  ruprcauntca  by  a  noiuul  pxouability  Uonslty  diatributicn 
defined  by  the  lolutlou  (^)  where  y  and  ai-e  the  die  an  and  standard 
deviation  and 


■4 


<0 


Invea tit^u tiuna  uf  conu  auiicut ion  pi  oblenxa  uaauciatod  with  random  noiae, 
which  has  ii/iny  uUvioua  aimllni  itiea  to  turbulence,  have  shown  that  normal 
dlatributiuJia  are  froquoutly  uiicountoied.  liice,  in  hla  'Math.  Anal,  of 
itandum  hoiau ”  haa  shown  that  for  a  lijtuar  ayst^m  the  shot  effect  in  a 
vacuum  tube  gives  rise  to  a  noise  current  which  haa  a  normal  diatiibution 
of  cuiront  intunaity.  Invoatigation  of  fluid  turbulence  frequently  yields 
normal  distributions  oi  velocity  f luctuationa . 


:jtutlaticully ,  whcnevoi  certain  conditions  are  mat,  the  distribution  in 
question  cun  be  consldei'ed  as  eighther  a  noriml  distribution,  of  aoioe 
skew  dlatributiun  which  can  bu  transfurntid  into  a  function  of  a  normal 
distribution.  V^hen  these  conditions  are  nut  standard  ^uethods  are  avail¬ 
able  for  toating  the  validity  of  all  stutintlcal  hypothesea  and  aaauniptiona . 


Meuaureiionts  indicate  tiiat  over  the  frequency  range  of  interest  the  power 
apectrum  of  a  tmospiiei  ir  turbulence  nuy  be  well  uppi’oximatud  by  simple 
analytic  expressions  such  as  have  been  used  in  wli*d  tupnel  atuaiea  of 
isotropic  turbulence.  Also,  in  the  theory  of  random  processes,  relations 
have  been  derived  between  peaK  counts  and  the  associated  power  spectra. 
These  I’elations  apply  to  the  case  oi'  a  stationax'y  uaussiau  oi’  hormal 
random  process,  where  the  term  uaussian  designates  a  process  characterized 
by  a  dausaian  distribution  for  the  amplitude  ol'  the  disturbances  us  wall 
as  for  its  time  derivatives.  Actual  xesulta  suppoi't  such  us  assumption. 


bincii  tho  Intonaity  of  turbulence  vai'loe  with  weather  conditions,  the 
operetionul  (^ust  history  is  considered  to  be  a  nonstetionery  Gaussian 
process  varyln^^  only  in  Intensity  or  rius  gust  velocity;  and  the  problem 
of  spooifyin^  the  gust  histoz’y  is  reduced  to  that  of  specifying  the  pro> 
bability  distribution  of  the  rms  gust  velocity. 


CUi-LIUTIVh,  H\Ub/a3lLrn  Ul:iThlbUriOf<  cy  lu-b  GLBT  VKLOCm 
^JUA  V:?.  LbAi-  DATA 


On  tliu  Dualu  01’  iticuu  cuiialduruiluiia .  utchniquoa  been  davelopod 

fur  the  eatl.uiCion  of  tiiu  probability  diutributions  of  the  rioa  accalaratlon 
and  tlto  rnLS  t/uai  value ity  from  data  on  \xuK  accelarations •  These  basic 
tachniquaa  dovalopod  by  Dr.  Harry  Proaa,  and  uaaociataa  at  NAHA  were  applied 
to  available  pouK  count  atatistica  from  b-47i  D-32,  and  KC-135  aircraft 
and  the  assoc iatod  probability  distributions  of  root  man  square  (just 
velocity  are  derived. 

Limitations  and  application  ol'  these  results  to  calculation  of 
t^ust-load  and  other  airplane  response  historiesi  i;:)ee  figure  1) 

1.  This  fcjraph  i-oprosunts  hours  of  b-q/,  b-52,  and  KL’-135 

flitjht  titna  abovo  lO.UUU  feet.  In  do terminiuti  the  transfer  functions, 
avorotio  values  of  airplane  and  operational  characteristics  Wei's  used. 

The  broken  curve  is  sui^orimposed  as  a  basis  of  comparison.  It  represents 
HAtA's  results  in  transport- typo  oporutions  from  jU.OOU  -  5D,0OU  feet. 

2.  The  cumulative  probability  distributions,  shown  wore  obtained 
by  inlctiratinti  the  probability  distributions,  ana  define  the  proportion 
of  total  flijilit  tiiio  spent  in  tui'bulencu  excoediiiji  given  values  of  (TtA.. 

J.  The  description  of  the  gust  experience  in  this  form  is  directly 
applicable  ^o  loud  calculations  for  other  airplanes  In  similar  opeiations 
by  reversing  the  proceduies  used  in  obtaining  tiiese  lesults.  however, 
direct  application  of  these  results  would  apply  oiily  to  similar  operations. 

In  order  to  obtain  results  that  are  mui'e  flexible  and  applicable  to  arbitrary 
flig^^t  plans,  it  would  be  desiruble  to  determine  tlio  vai'iations  in  these 
dis  tributionu  with  alt^uude,  wenthei’  condition,  and  perhaps  georgraphy. 

Perhaps  the  iiiosL  Important  points  to  be  noted  from  these  graphs  are 
the  relatively  lui'ge  amount  of  tiiri)  spent  in  essentially  sntooth  air  at 
the  iiigiier  altitudes,  and  the  relatively  large  ancunt  of  tiiic  in  ligiit 
to  severe  turbulence  at  tiie  lowest  altitude  biacive  t  (e.g.(^  is  greater 
than  2  feet  per  second  25/b  of  flight  time). 


Oivuji  Lilt)  V  iri  lilt  IL  l:i  iioL  yi;  L  poaiibla  to  duLuj'iailio  With  pi'wclaion 
tiiu  uiioi't;,  ui  LtiiLui  ..1.0  uii  aii-uciuica.  I!  il.a  vurlaliuna  aro  duo  to 
wuvu.i  witii  w.ivv)  luitr.Lho  01  1  mile,  the  uM'oclo  will  L>«  quite  dil'foreot 
than  il  iho  wovu  icii^tiio  wuro  100  fool,  ^uiuj  charac toriat lo  llRu  a 
‘acuio*  u!'  tuibuluiicu  io  iiuodod.  liuru  ,  vuriuiicu  ‘apoctra'  will  bu  uaod 
I'or  thia  piii  poou .  A,.,aiji,  tliouu  ohow  whul  rt'action  oi'  variauco  la  con- 
irlL/utud  by  vorioua  ii.lorvalu  of  I'ruquoncy  or  wave  loiif^th.  An  analO(^ 
la  111’)  oj  Ileal  ai^ietriLU  wMeli  aliuwa  that  I'l'uCtiun  ol'  the  atiort^  of  lit^it 
coiiua  f ru.u  vtiriouB  intorvula  of  fruijuojicy  or  wave  length.  Dtatiatlca 
tiuLliui-ud  lo  Jato  wolu  dui’ivud  from  vrih  rucorda.  i'ho  contribution  of 
utiiiodptu^ rie  lui'bulojieu  lo  ttiu  lullt/uo  ol'  aircraft  structurea  haa  not 
boon  coiuploloiy  •iiuilud  down"  to  dulo  -  in  any  rit^oroua  mathainaticul 
auliau.  lU>wiJVur,  aa  al.uwn  In  the  litarutura  and  in  lha  practical  uxparl- 
uuiit,  ti^od  atatiatlcbl  data  cojiibinud  with  Power  bpactrul  ilensity  Analysis 
Ima  proven  tjuitu  useful  in  dulineutinrf  re^ilona  of  lartjo  turbulence  Intensity 
which  do  cauao  coiia ide ruble  dama^je  to  airoruft  structuroa. 


